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Abstract 

A quantum field model of helical MHD stochastically forced by gaussian 
hydrodynamic, magnetic and mixed noices is investigated. These helical noises 
lead to an exponential increase of magnetic fluctuations in the large scale range. 
Instabilities, which are produced in this process, are eliminated by spontaneous 
symmetry breaking mechanism accompanied by creation of the homogeneous 
stationary magnetic field. 

1 Introduction 

Quantum field theory method including renormalization group (RG) approach has 
been successfully used for the theoretical explanation of various phenomena in de- 
veloped turbulence (see P and wherein references). 

In this paper the quantum field model of helical MHD is investigated. As a 
starting point we consider the Navier-Stokes equation for the velocity field and the 
equation for magnetic field which are driven by gaussian random forces with a given 
2x2 matrix D of the hydrodynamic, magnetic and mixed noise correlators, respec- 
tively. 

In ref. j2] (see also [HI) the multiplicative renormalizability of the quantum field 
model of non-helical MHD turbulence has been proved and the RG approach has been 
applied to study the asymptotic behavior of the model considered. The existence 
of two infrared-stable fixed points has been established. These points govern the 
two critical regimes: the magnetic regime and the kinetic one (the later being of the 
Kolmogorov type). 

The critical properties of the helical MHD are not known in the case of an arbi- 
trary noise matrix D. To provide the multiplicative renormalizability and consequent 
application of RG it is necessary to extend the theory adding the extra dissipative 
terms with new helical Prandtl numbers |3]. Therefore, also a critical behavior of 
the helical MHD is more complicated. A priori, the existence of the former stable 
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regime of the Kolmogorov type is not clear. In the following the existence of the 
critical regimes mentioned above for ordinary MHD is demonstrated for the helical 
one. 

There is an additional problem in the helical MHD: the instability of the theory 
which is induced by the exponential increase of the magnetic fluctuations in the 
large scales range (see 0, for example). The elimination of this instability leads to 
formation of a large-scale magnetic field known as the turbulent dynamo. Removal 
of the instability in quantum field formulation of helical MHD can be achieved by 
means of a nice and very well known spontaneous symmetry breaking mechanism 
followed by the creation of homogeneous stationary magnetic field. The special case, 
when only the hydrodynamic noise does not vanish, was analyzed in [H]. 

In this paper the value of spontaneous magnetic field c in one loop approximation 
for matrix D of noises in generic form. This value has been found from the conditions 
of overall exponential instabilities elimination in the steady state. 

The dynamo mechanism is accompanied by the appearance of an "exotic" term 
in the linearized equation for magnetic field. This term causes the linear growth in 
time t of the amplitudes of Alfven waves for small wave vectors k in the direction 
orthogonal to the plane of k and c. Due to the viscosity terms this growth is trans- 
formed into long-lived pulses of the type t exp(— i/Jt) exp(— at) with small a > and 
(3. 

2 The formulation of the problem 

The interaction of electrically neutral conductive turbulent incompressible fluid with 
the magnetic field with the magnetic field is described by the stochastically forced 
MHD equations "2\: 

dtv = z/Av - {vd)v + {hd)h -dp + F^ 

dth = z/'Ab - (v(9)b + (ba)v + , (1) 

together with the incompressibility conditions 

V ■ V = 0, V • = 0, V • = 0. (2) 

The first equation is the well-known Navier-Stokes equation for the divergence free 
velocity field v(a;) = {fj(x, t)} with the additional nonlinear contribution of the 
Lorentz force (p is a sum of both hydrodynamic and magnetic pressure per unit 
mass). The second equation for magnetic field b(a;) = {6j(x, t)} (it is connected 
with magnetic induction B by the relation b = 'B/^/AtFq, where g is the fluid density 
and n is its magnetic permeability) follows from the Maxwell equations for continuous 
medium. The magnetic diffusion coefficient u' is connected with the coefficient of 
molecular viscosity by relation z/' = uu with dimensionless magnetic Prandtl number 
(PN) u-\ 
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External random forces F'*', F'' are assumed to have a Gaussian distribution with 
< F >= and they are defined by 2 x 2 matrix of the noise correlators D =< FF >. 
The matrix elements are: the hydrodynamic D^^ noise, the magnetic D^^ one and 
the mixed D^^ one. 

The problem (P) is equivalent to a quantum theory with the doubled number of 
the fields $ = {v, b, v', b'} and the action functional [H [S]: 

, w'D^^'V h'D^^h' v'D^^b' b'D^^v' ^ 
S{^) = ^ + ^ + ^ + ^ + v'[-9tv + z/oAv 

-{yd)w + {hd)h] + h'\-dth + Moz/oAb - (v9)b + (b9)v] , (3) 

where v', b' are some auxiliary vector fields. Hereafter in the similar expressions, the 
integration over x, t and the traces over the vector indices are implied. As it is usual 
in QFT, the action © is considered to be unrenormalized with the bare constants 
marked by the subscript "0". The basic objects of the study are the Green functions 
of the fields $ (the correlation functions and response functions in the terminology 
of the original problem (J])). They can be determined as functional derivatives with 
respect to an external sources A = {A^, A^, A'^ , A^ } of the generating functional 
G{A) = J exp[S{^) + i.e., they are the functional averaged values of the 
corresponding number of the fields with a weight exp[S'(0)]. Here, D$ denotes 
the functional measure of the integration over the fields $ with all normalization 
coefficients. 

We have to choose a concrete form of D in the wave vector-frequency (k, u) 
representation. The noises are transversal for the incompressible fiuid. The action 
for helical MHD can possess scalar terms as well as pseudoscalar ones. Hence, the 
tensor structure of all noises is a linear combination of both tensor and pseudotensor. 
Then the correlators have the form: 

D): = Djl> = gy,k'-^'^-^^P%. (4) 

Here, P^^ = Pis + ipr£jsiki/k, where Pis = 6is — kikg/k stands for transverse projector 
and Eisi is Levi-Civita pseudotensor. Dimensionless real parameters p = {pi,p2,p3} 
satisfy the conditions | p |< 1, P3 <| P1P2 I- The scalar parts of the noises explicitly 
written in are in a standard power form [2]. The parameters go,gQ,gQ play the 
role of the bare coupling constants, and a, e are free parameters of the model. The 
value e = 2 corresponds to the Kolmogorov energy pumping from infrared region of 
the small k. 

3 Renormalization 

In a standard way, one can solve the primary infrared problem for the physical 
value e = 2 by the transfer to the region of small values e, where the ultravi- 
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Figure 1: One-loop Feynman diagrams which are UV- divergent. Only diagrams 
related to the Green functions < v'v >, and < b'b > (first and second line) can 
contain A — UV divergences 

olet (UV) divergences appear. They can be eliminated by the addition of the 
appropriate counterterms to the action Q 4J. The counterterms are formed of 
the superficial UV divergences, which are present in one-particle irreducible Green 
functions j^. The following 1-PI Green functions possess the UV divergences: 
< v'v >, < b'b >, < v'b >, < b'v >, < v'bb >. Corresponding diagrams are 
shown in Fig. 1 and Fig. 2 and counterterms have the form: i/v'Av, z/b'Ab, v'(b9)b, 
t/v'Ab and i/b'Av. The last two of them are not present in the primary action Q. 
For this reason, it is necessary to consider the extended theory with the additional 
cross dissipative terms viyv'Ah, wuh'Av and helical magnetic Prandtl numbers v~^, 
w"^. Besides UV divergences mentioned above which manifest themselves like the 
poles of e (e— UV divergences), another divergences proportional to the UV cutoff 
A can appear in the Green functions < v'v >, < v'b >, < b'v >, < b'b >. They 
acquire the form of Av' rot v, Av' rot b, Ab' rot v and Ab' rot b. These A — UV di- 
vergences generate the instability of the model, which causes exponential growth in 
time of the corresponding response functions. Therefore, their direct insertion into 
the action Q is not allowed and one has to find an effective way to eliminate them. 
One can make a natural assumption that finally the energy of the unstable large 
scale magnetic fluctuations must be transformed into the energy of large scale mean 
magnetic field. 

The model of the helical MHD under consideration describes steady state, there- 
fore it is reasonable to consider the new vacuum state with zero mean values of fields 
V, v', b', and non-vanishing time-independent mean field < b >= c 7^ 0. 

In a quantum field theory the appearance of non-zero vacuum value of field is 
associated with spontaneous symmetry breaking ^j, and as a standard, the value 
itself is determined from requirement of minimum of potential energy at the tree level. 
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Figure 2: Diagrams related to the Green function < v'bb >. 

In the case considered here the situation is more comphcated and rather technically 
different, because unstable A-terms appear at the next (one-loop) level, consequently, 
the nonvanishing value of mean magnetic field in the steady state must be calculated 
in this order of a perturbation scheme. By straightforward calculations and/or from 
the symmetry analysis of the given one-loop Feynman diagrams, one can find that 
only < b'b > contains A— terms. Here all A— divergences can be eliminated by 
means of the shift of b, namely b(x) — >■ b(x) + c, and, on the other hand, e— UV 
divergences are compensated by means of five independent renormalization constants 
Zj, i = 1, ...5 in extended model of helical MHD. As a result one obtains model with 
renormalized action: 

v'D^^v' h'D^^h' v'D^^b' b'^^^^v' 

\ 1 \ v'ldtv - 

2 2 2 2 ^ 

Ziz/Av - Z^vuAh + (v9)v - Z3(hd)h - Z3{cd)h] - (5) 
h'[dth - Z2Uiy Ah - Z^wuAv + {vd)h - (b9)v - (c9)v], 

where all parameters are renormalized couterparts of bare ones. The action (0) 
generates Green functions without divergences. In this case Feynman rules have the 
following form 
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^2 L+M+-S+V+' L+M+ - S+V+ ' 



'12 



(L+M+ - S+V+)(LM - SV) 



'12 



(L+M+ - S+V+){LM - SV) 

+ Df.ML+ + D^.VS+ - D\\VL+ 



where 



~ (L+M+ -,5+y+)(LM-,sy) 

(L+M+ -,5+y+)(LM-,5y) ' ^' 



L = —iuj + vk^ , M = — icj + vuk^ , 

V = vvk"^ — ^1 1 S = uwk'^ — i'j , 

J~)bv D^^ 

mn mn 1 

DZu = 52^'A;^-'-'^»(P^„ + ip2£mn;|), (7) 



with 



7 = c-k. (8) 

and vertices are defined by the expressions: 

v' ■ (v ■ (9)v = v[UjiVjVi/2 , 
v'-(b-9)b = viUjibjbi/2, 
b'-(b-a)v-b'-(va)b = ft^ii.zMi, (9) 

where 

tiji = + kAj) , % = i{kj5ii - kiSij) . (10) 

Using Feynman rules defined above one can immediately calculate the diagrams 
which are shown in Fig. 1 and Fig. 2. To extract the A— terms it is enough to keep 
only linear in wave vector part of the diagrams. The response functions < v'v >, 
< v'b >, < b'v > do not possess any linear divergent terms. A-divergent part and 
term ~ c of the response function < b'b > at the one-loop level are of the form: 

< b'Jjj > ~ ikm£imi{gpiCi + g'p2C2 + ^''PaCa) x 

.371 



X 



I'ASji — |c| 



'1 + ?iP — (v + wV 



8 \\ {l + uy{u-vw) ^ 



{6ji + ejei) 
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where e = c/|c|, Ci = {w{v + w) — u{l + u))/^, C2 = {1 + u — v{v + w))/^, 
C3 = 2{uv — w)/^ and ^ = 6(1 + u)^{u — vw)tt'^. From requirement of vanishing of 
A— term in (jlip one determines the value of spontaneous field 



Stt^ {\+uY{u — vw 

One can see from this equation that magnitude of spontaneous field is independent 
of coupling constants g,g',g", which characterize an intensity of random noises, and 
helical parameters pi, P2, P3- In such a way we obtain the renormalized Green func- 
tions which are finite as A — cx) formally, as it is usual in the field theory. But in 
real problems a natural maximal cutoff exists. In the developed turbulence the Kol- 
mogorov dissipative length = A~^ plays the role of a minimal scale. This length 
can be expressed in terms of basic phenomenological parameters - viscosity u and 
energy dissipation rate e; Id = u^^^€~^^^. Then from (fT^ one obtains |c| ~ {ueY^^ 
and it determines the order of magnitude of the spontaneous field c. 



4 Corrections to the Alfven waves 



In order to understand the role of the last term in (|TT|) let us consider the linearized 
MHD equations which follow from (0) in infrared limit of small wave vector k: 



dth 



—wvk"^ + Z7)v — uuk'^h + zx[k x e](b ■ e) 



(13) 



Here x = H9PiCi+g'p2C2+g"p3C3)\c\^{{l + u)^ - {v + wy)/{u - vw)/ {8-k{1+u)) 
and 7 = (kc). To solve Eq. p3|) we define a basis of orthonormal vectors n = k//c, 
1 = (e — n cos S) / sin 5, m = [n x e] sin 6, where 6 is the angle between vectors n 
and e. The transversal fields v, b can be decomposed with respect to the vectors 
V, b: V = vi\ + Vmm, b = 6^1 + bmm. Time dependent amplitudes vi, Vm, bi, bm satisfy 
equation (fT^ and are given by the following expressions: 

1 r 



vi 
bi 



cie^'\X2 -d) + C2e-'^*(Ai - d) 



cie 



(14) 
(15) 



X 



+ 



2Ai — a — d 
{e^^^ (^-cibe 



-C2be 
1 



2Ai 



2Ai ~ a — d 



2Ai - a 



e 



- + + C3 (Ai - a) - Cib^ + 

- + C3 (Ai -d)+ cS) ^ , 
n — fl \ 

L + 2t (Ai - d) 



2 \2\i-a-d 



(16) 
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- + - + H (2X7^ + 1 + 2t (A, - a)) 

+ C3C + C4(Ai -a)j^ , (17) 

where Ci, C2, C3 and C4 are constants of integration and we have defined 
1 



Ai,2 = -{a + d± ^{a + dy - A{ad-bc)) , 

a = —uk"^ , 

b = i'y — wk"^ , 

c = ^7 — wvk'^ , 

d = —uuk'^ , 

3 X 

e = -ik sin^ 6—— T\c\{gipiCi + g2P2C2 + gsPsCs): 



4(27r)^ ^1(1 + m) 



where 



'(1 + u)2 - (t; + ty)2 

M — VW 
2/ 



^1 = (1 + u) {u — i;w) . 

The functions Ai 2 are complex with negative real parts for the physical values of 
Prandtl numbers and they suppress the linear increase of amplitudes. 

Therefore, the last term in (jlip results in the appearance of specific long-lived pulses 
of Alfven waves. They are orthogonal polarized with respect to the spontaneous field 
and try to restore the isotropy broken by them. 

5 Short remarks to the RG analysis and critical 
regimes 

After successful elimination of all UV divergences we can use the renormalization 
group procedure and arrive to the set of RG Gell-Mann-Low equations for five in- 
variant charges: 

^~^ds^ = (^g^{9{s),e) , gi{s) \s=i= gi g^ = g,g\u,v,w, (19) 

where s = p/m (m is scale setting parameter). RG /5-functions are expressed via 
renormalization constants Zi, which have been calculated in one-loop approximation. 
Note that the /5-functions are finite in the limit of e — 0. 

The Gell-Mann-Low equations (|19p have been solved numerically for the various 
initial values of the invariant charges gi. It provides the possibility to analyze the 
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attracting regions of infrared fixed points. There arc two infrared-stable fixed points: 
the Gaussian (7* = and nontrivial g* ^ 0, u* = 1.393, g'* = u* = v* = w* = 0. 
The latter provides the existence of asymptotic critical regime of the Kolmogorov 
type. The RG approach improves expressions of simple perturbation theory and 
leads to the replacement of the original charges by the invariant ones, and, in critical 
regime by their values at the corresponding fixed points. Keeping in mind that c 
is connected with magnetic induction B by the relation c = B/^y4^^p|J (remember 
that p denotes density and /^-permeability of fluid) in above asymptotic region the 
magnitude of mean spontaneous magnetic induction is equal to Y\/^(^^y^^- 

6 Conclusion 

In the paper the correlation and response functions of velocity and magnetic fluc- 
tuations have been studied. Generally, these functions possess singularities, which 
can be eliminated by a proper renormalization procedure. As a result, RG equa- 
tions have been obtained and their solution have been found in the range of small 
wave numbers. This solution corresponds to the famous Kolmogorov scaling law. In 
helical MHD, where mirror symmetry of the system under consideration is stochas- 
tically broken, the non-vanishing mean magnetic field is spontaneously generated. 
This phenomena is accompanied by the linear growth in time t of the amplitudes 
of Alfven waves for small wave vectors. Due to the viscosity terms this growth is 
transformed into long-lived pulses of the type texp{—if3t) exp(— at) with small a > 
and p. 
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